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In electrical protection, there is a method of electrical protection of buildings 
against atmospheric discharges called the electro-geometric method or the 
rolling sphere method. So far, it is possible to achieve the implementation of 
this method graphically, that is, representing through plans and technical 
drawings, the protection conditions of the analyzed structure and obtaining 
from these graphic representations the protection parameters with the 
consequent errors caused by the scales and dimensions of the work plane. In 
the present work, a mathematical model is obtained that allows, using 
specific calculations, to analyze the dynamic behavior of a protection system 
against atmospheric discharges without worrying about the limitations given 
by the scales and planes. The set of equations obtained in the model allows 
us to determine the different parameters that define the protection system 
against atmospheric discharges (lightning) without depending on the 
graphical representation of the system's topology. 
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1. INTRODUCTION 

The electrical transmission, sub-transmission, and distribution systems must guarantee the safety 
and reliability of the electro-energy system in case of sudden changes in voltage and current, especially when 
there are atmospheric discharges. It is a common practice that in the protection systems against atmospheric 
discharges or lightning, Franklin-type lightning rods are installed in the corners of the buildings with their 
respective downspouts. However, on many occasions, there is no clarity of the area protected within the zone 
of lightning protection [1]-[3]. Various investigations have been carried out on the protection systems against 
atmospheric discharges and the different methods for their calculation and dimensioning. McDonagh and 
Klopotan [4] a 3D algorithm is designed to apply the rolling sphere method, analyzing the structures found 
within the area to be protected, as well as a combination of conical, spherical, and cylindrical sections. An 
analysis of a lightning protection system in conjunction with a grounding system is discussed in [5] to protect 
offshore oil platforms (OOP). 

A study that analyzes the overvoltages caused by atmospheric discharges in medium voltage lines is 
carried out in [6], determining the optimal location of the lightning rods to protect the transformers installed 
and in operation. Parise et al. [7], a study of a building in a complex of structures is presented, highlighting 
the worst scenarios of the impact of lightning in the implemented protection system. An investigation carried 
out in the pulp and paper industry is detailed in [8], where installing an external lightning protection system 
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using a single pole and the independent mast is described. The study details the best location of the lightning 
rod system, concluding that the extended mast terminal and the rod terminal are installed in the upper 
structures of the recovery boiler, the chemical room, and the ovens. 

Although there are standards and techniques to determine the protection systems against 
atmospheric discharges [9]-[16], this research work focuses its main objective on the mathematical 
development of the rolling sphere method, also known as the electro-geometric method, to determine the 
equations that they govern the behavior of the protection system and thus precisely determine the safe 
protection area under the electrodes or lightning rod. Usually, the electro-geometric method is used 
graphically [17], [18]. Still, this work proposes to find a mathematical model that consists of a set of 
equations that allow the calculation and design of protection systems against atmospheric discharges using 
only geometry analytical, trigonometry, and algebraic calculus by the rolling sphere method. In addition, 
different case studies are analyzed to validate the set of equations of the protection system. 

The rest of this document is organized as being as. Section 2 offers an analysis of the methods most 
used in the calculation and design of lightning protection systems used in buildings. In section 3, the bases 
are created to obtain the mathematical model of the shielding system using the electro-geometric method and 
the rolling sphere, evaluating the proposal in two case studies. Finally, in section 4, the conclusions are 
detailed according to the results obtained. 


2. METHOD 

The use of protection systems against atmospheric discharges and their different typologies are vital 
in interconnected or isolated electrical systems and buildings of all socio-economic sectors worldwide [19], 
[20]. Overvoltage is increasing in voltage that can appear on electrical distribution lines, and that can cause, 
among other things, damage to equipment connected to the network. Overvoltage is produced by direct and 
indirect lightning discharges, disconnection of inductive loads, network switching, and defects [2]. A study 
on the current and voltage levels of the protection system of a residential building and a nearby power line is 
analyzed in [21]; the authors model the system installed in the city of Huta Poreby in Poland, concluding that 
the value of the current between the low-voltage side and the 500-meter-long power line was only 0.02%, so 
the impact direct lightning on the power quality is relatively small. 

Naccarato et al. [22] the impact that a lightning strike produced in a residential building in the city 
of Sao Paulo, Brazil, in 2016 is analyzed. The authors conclude that the value of the current of the first peak 
saturated the Pearson sensor and that the systems ray localization (LLS) did not detect it, possibly because it 
had a slower pulse than the other spikes later than if the LLS detected them. Liu [23] designed software to 
determine the height of lightning rods with the help of the Borland C** Builder tool; with this program, it is 
possible to determine the height of one or more lightning rods using the polygonal method. For the 
verification, they used several examples, concluding that after comparing and verifying the results, the 
software is reliable and greatly reduces the calculation load during the design of lightning rods. 

A study on lightning protection systems in heritage monuments is presented in [24], where the 
modified rolling sphere method (MRSM) is used. In the research, the location of vulnerable points that can be 
hit by lightning is analyzed and, therefore, the area of the lightning rods is obtained according to the rolling 
sphere method (RSM) by which the impact distance is in as a function of the value of the electric current; 
however, the model ignores the dynamics of propagation and the geometry of the structure. Spunei ef al. [25] 
the damage that has occurred to the equipment and household appliances in an apartment building that has a 
protection system against atmospheric discharges is analyzed, concluding that it is necessary to update the 
standards of protection systems in buildings of more than 20 years old and placing protective screens near the 
electrical network to limit the induced voltage, in addition, they recommend installing the grounding 
conductor at a greater distance from the building in these types of constructions. 

Aslani et al. [26], by using the leader progression model (LPM) in 3D space together with an 
intelligent algorithm called teaching-learning based optimization (TLBO), propose a new method to identify 
vulnerable areas with the highest probability of lightning strikes. Using the TLBO algorithm, they define a 
new method to determine the critical current range and determine the most vulnerable points of the structure. 
The best lightning rod positions can be specified based on the simulation results. A sample of the asymmetric 
structure of height less than 60 meters is analyzed as a case study. The optimal locations and heights of the 
lightning rods are determined using the RSM and the proposed method, respectively. The simulation results 
validate the effectiveness of the proposed method, which makes it suitable for other complex structures. A 
study of the RSM and the Angle method [27], [28] to design protection systems against the direct incidence 
of atmospheric discharges on structures. This study shows, through mathematical modeling, that both 
methods are closely linked and that they can be used safely in the dimensioning of protection systems since 
they comply with the technical requirements demanded by the standards [29]-[33]. 
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3. RESULTS AND DISCUSSION 

Depending on the geographical location of the country in which the design of the protection system 
against atmospheric discharges is carried out, different ways can be found to carry out the design and 
calculation of the protection or shielding system [9], [10], but be it whichever method is used, in general, all 
are based on the distance traveled by the ray tracer to contact the surface to determine the protected area [11]. 
The distance traveled by the ray tracer is determined considering the risk level of the structure. It can be 
calculated using the standards International Electrotechnical Commission (IEC) 62305-2 [12], negative 
temperature coefficient (NTC) 4552-2 [13]-[16], or another international standard that contemplates it. 


3.1. Electro-geometric method 

According to the technical standard NTC 4552 [13] to analyze the action of atmospheric discharges, 
the electro-geometric model initially developed by [16], [17] is used and which has its application in the 
study of the shielding or protection provided by vertical rods and conductors horizontal. Considering the 
dielectric breakdown voltage between electrodes installed in large spaces, an empirical equation has been 
established that depends on the lightning return current (imax), proportional to the radius of a wait that rotates 
around the pick-up point, generating an area of protection (rsc). The equation that defines this radius is given 
in NTC 4552 by (1) [13]. 


Tee = Zimar + 30 (1 meee a) (1) 


Where: 
rsc: Radius of the rotating sphere around the pickup terminal and meters 
imax: Lightning return electric current and amps 

This expression can be simplified is being as (2). 


Tsc = 10 Gago” (2) 


Following the provisions of standard NTC 4552-2 [13]-[16], Table 1 shows the radius of the sphere of 
protection (rsc) according to the calculated risk level and the type of building where it is located installed 
protection system. In practice, to graphically determine the minimum height of the protection devices, 
circumferential arcs with radius equal to the impact distance rsc are drawn between the objects to be protected 
and the interception devices, for example, Franklin-type devices, in such a way that is tangent to the ground 
and to the targets that make up the arc. Any object in the line that describes the arc will be exposed to direct 
electrical discharges [18]. Figure 1 shows the concept of the electro-geometric model. 


Table 1. Risk levels according to the characteristics of the building 
Risk level rse (meters) Fse (feet) 


Level I 35 115 
Level II 40 131 
Level III 50 164 
Level IV 55 180 


Figure 1. Schematic of the electro-geometric model 


3.2. Rolling sphere method 
The RSM is a corollary of the electro-geometric method and consists of assuming a sphere of radius 
equal to the impact distance that rolls over the structures to be protected [21]. The radius of the rotating 
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sphere rsc is determined in the same way as in the electro-geometric method, using (1) and (2). However, it is 
widespread to take a radius of 150 feet or 45.72 meters as a reference. Figure 2 shows that the central or 
maximum circumference of a sphere has the largest radius, while any circumference parallel to it will have a 
smaller radius. When a sphere lies on two specific points that belong to the same plane that contains the 
sphere's center, it can be said that these two points belong to the central circumference. If these two points are 
not in the plane containing the sphere's center, the generated circumference will have a radius smaller than 
the maximum. When it is known, with the maximum diameter of the sphere, it is possible to simplify the 
mathematical model by posing the equations that can be used to determine the optimal distances between two 
consecutive points (electrodes of the protection system). 

Figure 3 describes the sphere's movement on the structure of a multi-level building; said sphere 
must roll contacting the tips of the pickup terminals and finally reach the ground without contacting the 
structure, thus creating what is known as a protected area. The rolling sphere method should not only be 
applied to the building but also across the width of the structure to be protected [22], [23]. Considering the 
diameter of the rolling sphere, the calculations can be simplified since when it contacts any of the pick-up 
points, it will do so at some point on the circumference. This simplification is valid because the sphere from 
its center to any point on its surface will have a radius equal to rsc. 


Sphere surface 


-X 
Ma Y 
E a 
Center circumference 
. Yo 
SX 
Protected zone 
Figure 2. Cross-section of a sphere that Figure 3. Schematic of the RSM 


describes a circle 


3.2.1. Sphere that travels between the ground and a collection tip 

When analyzing this case study, it can be observed that there are three critical physical points, the 
center of the sphere, the point where the sphere touches the ground, and the point where the sphere contacts 
the lightning rod, see Figure 4. By locating the coordinate axis in a convenient position for analytical 
calculation, the center of sphere C is identified as having coordinates C = (0, rsc). Also, the point where the 
rolling sphere touches the earth has coordinates D = (0, 0) and is the point where the coordinate axis of the 
system is placed. Finally, point P, which is located at the tip of the lightning terminal (lightning rod) with 
coordinates P = (Li, [hit+h2]), where L is the separation distance between point D (the point where the sphere 
touches the land) and the building, in meters; h; is the height of the building, in meters and hz is the height of 
the collector point, also in meters. These points are shown more clearly in Figure 5. 

The two-dimensional equation of the sphere is given by (3). 


Cay +Q-b? ax (3) 


Where: 
a: Value of the coordinate on the x-axis 
b: Value of the y-axis coordinate 
r: Magnitude of r;- according to the level of risk calculated 
By replacing point C, the center of the circumference, (3) becomes: 


(x = 0)? + = te)? = ne’, X? + (YY = Tre)? = Toe” (4) 
In (4) is used to determine the safety zones and values of the protection system. 


As this equation corresponds to the mathematical expression, it is possible to determine any 
coordinate (x, y) that belongs to said circle. 
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Solving for (4) for the x coordinate axis, the expression (5) is obtained: 
xX = Ņ 2Yfge = y? (5) 


Solving for (5) for the y coordinate axis and obtaining: 
X =V he? — X? + Te (6) 


Ground 


Ground 


Figure 4. Case study 1: ground and catchment Figure 5. Case 1. Coordinates of the analysis points 
terminal 


With (5) and (6), it is possible to calculate the points (x, y) that belong to the circumference and with 
which the safety parameters or the minimum distances that must be considered in the protection zone can be 
determined of the building. Point D shown in Figure 6 is the ground closest to the building that is most likely 
to be struck by lightning. Using (4), it is possible to determine its distance. 

If the coordinates of point P are replaced in (4), the expression that allows calculating the distance 
between the structure and the most probable point of impact is obtained, thus obtaining (7). 


Lm = J2(hy + hy) fse — (hy + h3)? (7) 


Where: 
Lm: Distance from the building to the most probable point of impact on the ground, meters 
hy: Height of the building to be protected, meters 
hz: Height of the selected lightning rod, meters 

After determining the maximum impact distance, the value of the physical design separation 
between the lightning rods of the protection system and the point of impact closest to the building must be 
calculated. For the design separation (La) calculation, a new factor is introduced, called the safety factor (fs), 
which determines the percentage that the design separation La can be increased before reaching the impact 
separation Lm. To determine La, (8) is used: 


La + fs) = Lm (8) 


Where, 0 <f < 1 
From the previous equation, it is possible to solve for the design separation La, and we obtain (9). 


La = 2 (9) 


~ fs) 


Another aspect to take into account is the protected separation (Ls), which is the difference between the 
impact separation value (Lm) and the design separation point (La), and can be calculated by (10). 


Ls = Lm — La (10) 
We also have to analyze the protected height (h,), the value of the existing separation between the 


ground, and the point of the rolling sphere where the earth and the lightning rod (lightning rod) make contact. 
In Figure 7, the height h, is represented, and it is also observed that the value of the x-axis = Lg. 
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ec 


Figure 6. Case 1. Coordinates of any point (x, y) Figure 7. Protected height hp 


To calculate the protected height hp, the value of the coordinate x is replaced in (5), obtaining the 
expression (11). 


hp EN Tsc? = (La)? + Tsc (11) 


Where: rsc is the sphere's radius determined with the calculation of the risk level. 

With the value obtained from Lm in (7) and using (4) and (5), it is possible to determine any 
coordinate of the circumference shown in Figure 6 to analyze whether a specific area or height is protected. 
Figure 8 shows the importance of the security zone covered by the shielding, where a structure located next 
to the building can be seen; this structure represents any element, for example, a water or fuel storage tank, 
an antenna, and a silo. As shown in Figure 8, point E is exposed to the direct impact of atmospheric 
discharges, so it is necessary to expand the protection zone to completely cover the building and the structure 
next to it. The expansion of this area is achieved by increasing the height of the lightning rod, which causes 
the sphere to recede away from the building and expanding the protected area, thus covering the structure to 
be protected. In Figure 9, it is observed how as the height of the lightning rod increases, the sphere moves, 
and the location of point E changes, contacting the sphere and leaving the object within the safe zone. 


ec i 


' 
I 
I 
I 
I 


Figure 8. Unprotected structure Figure 9. Protected structure 


3.2.2. Sphere supported on two points located at the same height 

The second case study is represented in Figure 10 and consists of two lightning rods separated by a 
defined distance. These lightning rods have the same height and contact points A and B with a sphere with 
center C and radius r;.. Points A, B, and C can be analyzed as an isosceles triangle, as shown in Figure 11. 
The two lightning rods have a height of hl and are separated by the distance L1. The sphere's point that 
contacts the plane where the lightning rods are located is at h3. The distance between points A and C and B 
and C is rsc, and the perpendicular distance between point C and line segment AB is h2. 
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c (0, [h1 + hA) 


AC-L1/2, hT) B (L172. h1) 


x 


pu — i = 
}-—— 1/2 —»j +L. 12 


Figure 10. Sphere supported on two lightning rods Figure 11. Case study 2. Sphere contact zone 
that are at the same height 


In the system shown in Figure 11, the coordinates of the different contact points can be observed; 
point C, which is the center of the circle, has coordinates (0, h)+hz). At point A, which is the tip of one of the 
two lightning rods, this coordinate is given by (-L//2, hı), and at point B, it is given by (L,/2, hı). The value of 
h2 can be calculated with the right triangle formed between the midpoint between A and B, the center C, and 
point A. Applying the Pythagorean theorem for this triangle, we obtain (12): 


re? = (2) + hy? (12) 


Solving for h2 from the previous equation gives (13): 


hy = ree? (4) (13) 


Working on the two-dimensional equation of the sphere shown in (4), we have: 
(x -—a)?+(y—b)? =r? (14) 


Where a is the value of the x coordinate, b is the value of the y coordinate, and r is the magnitude of rse 
depending on the level of risk calculated. 
By replacing point C, which is the center of the circumference, (1) becomes: 


(x — 0)? + (y — [hy + h2])? = i Fag O - hee h2])* z Tee (15) 


In (15) is used to determine the zones and safety values of the protection or shielding system. This 
expression makes it possible to decide on any value corresponding to the coordinate axis (x, y). 
From (15), the coordinate x can be solved, obtaining: 


x = Jre? — O — [h + hp)? (16) 


To obtain the y coordinate, we have: 


Y ZNT? =x" + (hy + h2) (17) 


These equations make it possible to obtain the coordinates of any specific point on the 
circumference. To determine the value of the separation between the structure and the rolling sphere, 
represented by h3 in Figure 11, the general expression given in (18) can be used. 


(hy + h2) = fsc + hz (18) 


From the analytical point of view, it is also possible to determine the point h3, since when the 
coordinate x = 0, then y = hz. 
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By replacing x = 0 in (17) we obtain that: 


y= af Gee + (hy + hz) (19) 
So that: 
Y = Enc + (hi + hz), hg = Tsc + (hy + h2) (20) 


In this case, two solutions are obtained for x = 0, one for +r;., which would be the upper point of the 
circumference, and another for -rsc, which would be the lower point, that is, A3, so the final equation is given 
by (21). 


hz = — Ts + (hy + h2) (21) 


It should be noted that (18) is the same as (21) obtained by two different methods. Substituting (13) in (21) 


1S: 
2 hh)? 
hy + Jre? - (4) — tee = hs (22) 


In (22) allows finding the separation (h3) between the rolling sphere and the structure that is being 
protected only with the basic design parameters, which are the height of the lightning rod (/,), the radius of 
the protection sphere (rsc), and the separation between the lightning rods (L,). An essential aspect of this case 
study is to determine the circumference equation when the maximum separation between lightning rods 
occurs, that is, to determine the maximum separation distance between lightning rods in which the rolling 
wait would contact the protected structure. As shown in Figure 12, the center of the circumference of the 
rolling sphere is located at the coordinate C = (0, r;-). By replacing this center in (1), which is the general 
equation of a circle, we get (4). 


Lm 


D EB EB Ga ba A Fa 


Figure 12. Maximum separation of lightning rods in a protection area 


4. CONCLUSION 

With the mathematical development carried out, it has been possible to create a new mathematical 
model that offers the possibility of calculating the different parameters of the lightning protection system 
without depending on the graphical representation of the topology of the system, which allows reducing the 
derived errors of the limitations of the scale on which the structure to be protected is graphed. This 
mathematical model is a useful tool in analyzing the different configuration cases that can occur in buildings 
that must be protected against atmospheric discharges without leaving unprotected areas. In the same way, it 
allows the designer of the protection system to have a precise dimensioning of the location of the different 
elements of the system. This mathematical model provides a tool that allows the design of computer 
applications (software) to calculate the protection system against atmospheric discharges using the rolling 
sphere method and considering the safety criteria according to the current technical regulations of each 
country. In addition, the change in the design variables of the system allows modeling different 
configurations to determine the most suitable and efficient one. 
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